Resonance phenomena of a soliton-like extended object in a bistable potential 
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We investigate the dynamics of a soliton that behaves as an extended particle. The soliton 
motion in an effective bistable potential can be chaotic in a similar way as the Duffing oscillator. 
We generalize the concept of geometrical resonance to spatiotemporal systems and apply it to design 
a nonfeedback mechanism of chaos control using localized perturbations. We show the existence of 
solitonic stochastic resonance. 
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The soliton dynamics in inhomogeneous media fl) - fl, the spatiotemporal chaos (J], Q- Jl0| , the chaos control |pd|| - 
fl5f and diverse types of resonance phenomena |]l6|- |)22f have become the object of intensive study in recent years 
(although rarely these phenomena have been studied simultaneously). 

In the present letter we investigate a solitonic model perturbed by inhomogeneous forces jf|- |J, for which the exact 
solution describing the stationary equilibrium solitonic state can be obtained. Furthermore, the stability problem for 
this solution can be solved exactly. The effective potential for the soliton motion can be bistable depending on the 
. system parameters. Using an additional external periodic force we can have the soliton jumping chaotically between 
the two potential wells (like in the Duffing oscillator) p3| . We prove that changes of the soliton shape and the 
wave-form of the external perturbations are very relevant to the soliton's dynamics as a whole. We generalize the 
concept of geometrical resonance p2| to spatiotemporal systems. This concept can be used to design a nonfeedback 
mechanism of chaos control 113]- DM. In our case, this mechanism can be applied in a localized way in space. Finally, 



we show the existence of solitonic stochastic resonance (SSR). 
I/"") . We are interested in equations of the type: 
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OS] 4>tt - <t>xx + R((f>,<t>t,x) - G{4>) 



F(x)+P(<l>,<k)q{x,t), (1) 



where G(<f) = — dU g^ ; U{4>) is a potential that possesses at least two minima |24), [p5| ; R(<j>, <f>t, %) is a dissipative 
+1j . term, which in general can be nonlinear Q and the damping coefficient can depend on x; F(x) is an inhomogeneous 
force, and P(<p,(f>t)q(x,i) is a general temporal perturbation. Many systems are described by this type of equations, 
O |- including charge density waves, Josephson junctions, and structural phase transitions j2f|. 

Problems related to Eq. (^) present extraordinary mathematical difficulties. Even if we are to solve (for case 
•rH . P{<j),(j)t)q{x^t) = 0) the stability problem of a stationary soliton placed on an equilibrium position created by the 
inhomogeneities |^], this task requires the making of big approximations The most common approximation is 
to consider the soliton as structureless point-like particle J| (in this case the coordinate of the soliton center-of-mass 
and its velocity are the only dynamical variables). Although this is a valid approximation in many cases, our work 
|0|- J5| has shown that, in general, the internal dynamics plays a fundamental role. 

Besides the applications of the 4 -equation in phase transition theory |26|] , it is a model system for topological 
solitons in general. Let us consider as an example the perturbed </> 4 -model: 

<t>tt-<Pxx+1<f>t- \^-<j?) = F{x), (2) 



where 
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F(x) = F^x) = ]^A{A 2 - 1) tanh(Ba;) 



1 . , „q .os smh(Bx) 
-A(AB - A 2 ) 



2 cosh. 6 (Bx) 

The inhomogeneities are chosen in order to have some important properties. The exact stationary solution for the 
soliton equilibrated by the inhomogeneities in the point x = can be obtained <pk = Atwh{Bx). The stability problem 
of this solution can be solved exactly. The force F\ (x) possesses the interesting property to be topologically equivalent 
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(in the sense of Catastrophe Theory p^ ) to a pitchfork bifurcation, allowing us to have an effective potential (for 
the soliton) with one or two wells, depending on the system parameters. This system is generic and structurally 
stable. Therefore the results obtained in this letter can be generalized qualitatively to other systems topologically 
equivalent to that described by Eq.^. The force F(x) permits to observe important phenomena which occur with the 
participation of the internal soliton dynamics including the appearance of a great number of internal modes and the 
soliton destruction due to the instability of the shape mode in some special situations 0-0]. 

The stability analysis 0- || which considers small perturbations around the soliton (4>(x,t) — 4>k(x) + f(x)e xt ) 

leads to the eigenvalue problem Lf = Tf, where L = —d 2 + {^A 2 — \— 2 cotS(Bx) ) > ^ = ~(^ 2 + T^)' 

The eigenvalues of the discrete spectrum are given by the expression T n = — ^ + B 2 (A + 2 An — n 2 ), where A(A + 1) = 
The stability condition of the equilibrium point x — (for the soliton) is defined by the eigenvalue of the 
translational mode (fa(x) = cosh~ A (_Ba;)): B 2 A — i > 0. A global topological analysis and the investigation of 
the infinity j^] allow us to have a complete qualitative information of the soliton dynamics. When the equilibrium 
position for the soliton is unstable and the absolute value of the eigenvalue corresponding to the translational mode is 
sufficiently high, the first shape mode (fi(x) = sinh(Ba;) cosh~ A (i?x)), and even other internal modes, can be unstable 
too, producing the soliton destruction p|. We have at least two interesting situations: If A 2 > 1 and 2AB 2 > 1, then, 
there exists only one stable equilibrium point for the soliton (xqq = 0). If A 2 > 1 and 2AB 2 < 1, a pitchfork- like 
bifurcation occurs and we have two stable equilibrium points (xoi < 0, xq 2 > 0) (the point xoo = becomes unstable). 
It is important to stress that the bifurcation does not occur at the point we would expect when considering the soliton 
as a structureless point particle (in this case the number of equilibrium positions depend on the number of zeroes of 
force F(x)). 

If the system, in the bistable case, is perturbed by an additional periodic force fitted to the shape of the translational 
mode {§, §] (in (Eq||) we put P{(j>,<t>t) = 1), then we can have similar situations to those of the Duffing oscillator 
P3fl . But remember that in this case what is jumping between the potential wells is an extended object with a very 
complicated internal dynamics. On the other hand, for sufficiently high values of |To| for the unstable position, the 
soliton can bifurcate § in an antisoliton (which would feel the middle position as a stable one), and two solitons, 
each of which would move towards one of the wells. 

The recently formulated concept of geometrical resonance [ p2[ (GR) can be very useful in this context. Generalizing 
this concept for spatiotemporal systems we define (f>Gn(x, t) as a GR solution of equation (Q) if it satisfies the condition 

R{(j)GR, 9 t (f>GR, X) = P(<pGR, 9 t (f>GR)qGR(x, (3) 

In this case there exists (local) conservation of energy: 



H 



1 ( d(j)GR\ 2 , 1 / d(j)GR^ ' 



2 V dt J 2 V dx 

U{(j> GR ) - F(x)<t> G R + C 



dx (4) 



If we want to observe a GR situation, we should use a qGR.(x, t) such that it satisfies the GR condition (|^). 
Consider the Eq. (|l|) with R = yfo; P = 1; G{4>) = \(<j> - </> 3 ) and F(x) as defined in (|). 

We assume that the stability condition B 2 A — ^ > for the equilibrium position x — 0, holds. In this case we can 
write an approximate solution for (j)(x,t) (q = 0): 

^,t)=Atanh(^) + feooCOS ^V 0) ' ( 5 ) 

cosh (Bx) 

where Wq =Tq. 

Thus, for small values of hoo the perturbation 

u h 0Q j sin(uj t + 6 ) 

= ^Khb^) (6) 

satisfies the GR condition. This explains a whole series of experiments preformed in [pj, j^j. There the </> 4 -kink confined 
in a potential well created by the inhomogeneity F(x) was perturbed by the time-dependent force (^|). The authors 
have got resonances of the translational mode practically without deformation of the kink profile. Chaotic behavior 
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is obtained only for high values of the perturbation amplitude, when (||) is not a GR solution anymore. Meanwhile, 
if we use a different time-dependent force (e.g., q(x,t) = pocosut) which does not satisfy the GR condition, we will 
obtain irregular behavior in time and space with much smaller amplitudes. 

The GR [^2| provides a mechanism for nonfeedback control of chaos. 

Suppose we have the following equation 



F(x) 



cos(u)t) 



We assume that for certain initial conditions the Eq. (|7|) is in a chaotic regime provided F c {x,t) = 0. 

We are left with the problem of eliminating the chaotic motion using control F c (x, t). Additionally, it is expected 
that the control driving term is small and localized in space. In order to do this we should select the control term in 
such a way that the solution will be sufficiently close to a T' periodic GR solution. 

We can choose the control function in the form 



F c {x,t) 



Using the concept of "almost adiabatic invariant" 



g c cos(w c i 



cosh A (Ba;) 



(8) 



|22| , we arrive at a condition for the control parameters 
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dt 

b GR cos(w c i 



Pq- 



•>GR 



dt 



cos(ui) 



+ C ) 



dt cosh 2 (Sir) 



d.r 



(9) 



The fact that we have taken (||) associated with the translational mode (which approximately satisfies the GR 
condition) always allows us to find control parameters (with small g c ) in order to suppress chaos. An Arnold-like 
tongue structure, similar to that observed in |||] can be obtained here. In Fig. 1 we show the results of the application 
of the localized nonfeedback control. 

Note that the validity of the concept of GR as a mechanism for chaos control is not limited to the case in which 
the unperturbed equation can be solved exactly. Using some topological analysis it is possible to guess the "shape" 
of the GR solution and then, through the concept of "almost adiabatic invariant" ^j"^ \ ~ O^one can obtain 

the conditions such that the solution will be close to a GR solution inside of a "mode-locked" tongue. 

The system we have presented in this letter with a soliton-like extended object moving in an effective bistable 
potential can be also useful for studying other phenomena, e.g. the spatiotemporal stochastic resonance (SR) for the 
motion of an extended state with complicated internal dynamics. 

Consider the equation 



- <t>xx + l<f>t - \{<P - 3 ) = F{x) + q(x, t), 



(10) 



where 



F(x) 



B(AB 2 - 1) tanh[B(a; + x Q )} for x < 0, 







for x = 0, 



B{AB 2 - 1) tanh[B(x - x )} for x > 0; 



(11) 



[P cos(a;£) + i](x, t)] cosh _2 [B(a; + x Q )] 
for x < 0, 

[Po cos(wt) + r](x, t)] cosh~ 2 [B(x — x )] 
for x > 0; 



here 4B 2 > 1. The forces F(x) and q(x,t) are defined by Eq. ( fLlpi! ) in such a way that there are two equilibrium 
points for the soliton and the motion in each well is very close to the GR condition. In the absence of the white noise 
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r](x,t) (< f](x,t) >= 0, < rj(x,t)ri(x' ,t') >= 2DS(x — x')S(t — t')), the periodic force above is unable to make the 
soliton to jump between the wells. 

When we switch on the noise, it is possible to observe a maximum in the graph of the Signal-to-Noise-Ratio 
(SNR) versus D (see Fig. 2). Here as the "signal" we take the time series of the soliton center-of-mass coordinate: 

J x<p^.dx 

x c . m , — ~l/2 • The SNR is measured following the traditional method 113] . 

J -1/2 x 

The maximum synchronization is obtained with a signal associated with the translational mode. These results allow 
to predict the optimum entertainment of the noise by means of a periodic signal, not necessarily a simple sinusoidal 
signal but a more complex spatiotemporal function. The stochastic resonance will depend on the characteristic shape 
of the kink and the wave- form of the perturbation (in time and space). A manifestation of this fact is the following 
phenomenon: if instead of the translational mode in perturbation (12) we use the first shape mode, then we do not 
observe SSR. 

Even when F(x) has three zeros, if we are in the presence of an extended object, the bistability is not a sufficient 
condition for the stochastic resonance. The extended object should "feel" the bistability and the internal modes 
should be stable in the vecinity of all the equilibrium points, including the central equilibrium point which is unstable 
for the translational mode. On the other hand, if the eigenvalue |r | corresponding to the unstable equilibrium 
position in the bistable potential is very high, then the first shape mode can be unstable too and the soliton will 
bifurcate in two solitons and one antisoliton. This is what occurs when the parameters in Eq. (^) are such that 
B 2 (4B 2 - 1) tanh(Bxo) > §§. The center of mass of this "three-particle structure" is always oscillating around x = 
and therefore, there is no SSR. 

The spatiotemporal stochastic resonance has been considered in some recent papers jl9), |2l}| . In particular, Ref. 
|jl9f deals with the 4 -equation. However, the SSR presented here is a completely new phenomenon: there the 
behavior of the field 4> is taken for the signal and the important bistability is in the potential U(4>), whereas in our 
work the relevant signal is the time series for the coordinate of the soliton's center of mass. In Ref. there is no 
inhomogeneous force F(x). The spatiotemporal stochastic resonance studied there does not depend on F(x). In our 
case the potential U(4>) is important for the existence of the soliton solution |24j] , but it is the bistability in x, created 
by F(x) which is the key for the SSR. And, as we have seen not for every F{x) the SSR exists. 

In general, we prove that changes of the soliton shape are very relevant to its dynamics as a whole. Additionally, 
the wave-form of the perturbation is crucial for all the resonance phenomena including the nonfeedback mechanism 
of the chaos control and the stochastic resonance. We believe that these phenomena are relevant to other systems 
where there is an extended state with a complicated internal dynamics moving in a nonlinear potential force. 
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FIG. 1. Phase portraits for the motion of the soliton center of mass in eq. (Q), F(x) is defined as in eq. (^) (B — 1/2). a) 
Chaotic motion in the absence of control (F c (x,t) = 0), 7 = 0.1, A = 1.5, Po — 1 and ui — 0.65. b) Periodic motion as a result 
of the application of nonfeedback control, ui c — 0.065, g c = 0.35 and 9 C = 0. 

FIG. 2. Two aspects of the solitonic stochastic resonance: a) There is a maximum of the SNR at a critical value of noise 
intensity, b) Sinchronization of the stochastic jumps with the periodic perturbation linked with translational mode of the 
soliton.(P = 0.04, u = 0.01, x = 2.5 and B = 0.7). 



FIG. 3. Spatiotemporal evolution of the kink-soliton at the stochastic resonance. Note that there is a minimal deformation 
of the kink profile 
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